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Abstract. We investigate the possibility of observing monochromatic neutrino lines originating from an¬ 
nihilation of dark matter. We analyse several astrophysical sources with overdensities of dark matter that 
can amplify the signal. As a case study, we consider mixed left and right handed sneutrino dark matter. We 
demonstrate that in the physically viable region of the model, one can obtain a prominent monochromatic 
neutrino line. We propose a search strategy to observe these neutrino lines in future generations of neutrino 
telescopes that is especially sensitive to dwarf spheroidal galaxies. We demonstrate that the presence of 
massive black holes in the cores of dwarfs as well as of more massive galaxies substantially boosts any 
putative signal. In particular, dark matter in dwarf galaxies spiked by IMBH provides a powerful means 
of probing low annihilation cross-sections well below 10~^®cm^s“^ that are otherwise inaccessible by any 
future direct detection or collider experiment. 


1 Introduction 

Indirect detection of dark matter (DM) is an important probe 
of the nature of DM. It relies on the search for annihilation 
products of DM particles, which can include photons, anti¬ 
matter and neutrinos. The usual strategy is to observe a sec¬ 
ondary flux arising from DM annihilations, for instance at the 
center of the galaxy or in point sources. Despite the astrophys¬ 
ical uncertainties involved, limits from e.g. space and ground- 
based Gherenkov telescopes set important constraints on the 
thermal cross-section of the DM annihilations into gamma rays. 
In terms of the neutrino secondary flux, high energy neutrinos 
produced by the capture and annihilation of weakly interact¬ 
ing massive particles (WIMPs) in stellar cores are one of the 
smoking guns for DM detection and are actively being 

searched for by neutrino telescopes. Several studies, e.g. [^|^, 
rely on this signature and help to determine future prospects 
for neutrino detection in supersymmetric (SUSY) models. 

Another search strategy for DM in astrophysical environ¬ 
ments relies on spectral features, which arise in some DM sce¬ 
narios, e.g. monochromatic lines that target the DM mass and 
serve as a smoking gun signature for WIMPs [^. Gamma-ray 
lines have been studied extensively (see this review for further 
references [^), while much less attention has been given to 
monochromatic neutrino lines [7||^, which however may have 
rather interesting features. 

In this paper, we study the prospects for detection of mono¬ 
chromatic neutrino lines from WIMP annihilations arising at 
tree level in certain DM scenarios. Neutrinos can be a partic¬ 
ularly clean signature of DM annihilations as they propagate 
freely and undisturbed to the Earth, while the other parti¬ 
cles undergo absorption and scattering processes before reach¬ 


ing the Earth. Indirect searches for DM are particularly sen¬ 
sitive to the astrophysical environment. This introduces sev¬ 
eral uncertainties in the limits derived using indirect detec¬ 
tion techniques. However, we take advantage of this fact and 
assess the impact of enhancement in the indirect detection 
signatures due to DM overdensities in diverse astrophysical 
sources. Among alternative astrophysical environments with 
an enhanced DM density, we consider galaxies hosting a su- 
permassive black hole (SMBH) and dwarf spheroidal galaxies 
(dSphs). We propose these astrophysical targets as a future 
detection strategy for DM particles at neutrino telescopes. We 
demonstrate that these searches can probe values of the anni¬ 
hilation cross-section well below 10“^®cm^s“^. 

The presence of a SMBH might have dramatic consequences 
for DM annihilation. Such black holes (BHs) may be common 
signatures of galaxy formation, at least as envisaged in many 
scenarios of structure formation [^[^, and are capable of 
boosting the DM density within its sphere of gravitational 
influen ce [12] . In particular BHs can act as a particle accel¬ 
erator [13], especially in the case of Kerr BHs. Eor WIMPs 
close to the BH horizon, the scattering may occur at almost 
inflnite center-of-mass energies [14}{16] , for example boosting 
cross-sections which are p-wave and usually suppressed at the 
present epoch [^, or exciting resonances and producing a 
forest of gamma-ray lines [18[|19] . In our case, the dominant 
cross-section for producing neutrino lines however is velocity- 
independent and we will not investigate velocity-dependent en¬ 
hancements here. 

We focus here on the more generic scenario in which the 
SMBH may increase the DM density distribution towards the 
center of the galaxy and form a so-called DM spike [^. In 
the standard picture, the DM density profile follows a power 
law, p oc r~^ in the inner region, with r the distance from 
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the galactic center and 7 ~ 1 the slope for a Navarro-Frenk- 
White (NFW) profile [^. However, very close to the BH, the 
DM density profile may rise very steeply. Assuming that the 
BH grew adiabatically, 7 is expected to instead lie between 
2.25 and 2.5, leading typically to 7 = 7sp = 7/3 in the very 
inner region. Such an enhancement of the DM density profile 
is referred to as a spike and enhances the luminosity of the flux 
originating from annihilating WIMPs. 


The existence of such a DM spike in the vicinity of SMBHs 
is still under debate. If the BH growth was not adiabatic, 
then the i nner D M density profile would behave instead as 
Besides, if the DM halo itself underwent a 


p (X r 
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merger, or if the BH did not grow exactly at the center of the 
DM halo, the inner DM halo profile would follow p oc 
thus considerably reducing the flux expected from WIMPs. 
Even if a spike could form with 7 sp ~ 7/3, the process of 
dynamical relaxation by WIMP scattering off stars could par¬ 
tially smooth the spike and lead to a much shallower slope, 
p oc , which would coincidentally correspond to a Moore 
. Constraints from a DM spike in terms of annihilat- 
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profile 

ing DM producing an electromagnetic flux have been discussed 
in [^ 26 . Use of a secondary neutrino flux to probe the DM 
spike has been initially proposed in [^, however monochro¬ 
matic neutrino lines have not yet been discussed in this con¬ 
text. 


Another interesting target for indirect detection of DM is 
dwarf spheroidal galaxies (dSphs) in the Milky Way Galaxy 
(MWG) [^, along with the Galactic Center. What makes the 
former really interesting is their low astrophysical background, 
high mass-to-light ratios and proximity to the Earth (see for 
example 28 ). To further enhance signals from WIMP anni¬ 


hilations, we study the impact of a DM minispike due to an 
intermediate massive black hole (IMBH) hosted in the inner 
region of nearby dSphs [^. Although the presence of the DM 
minispike is also a matter of debate and the answer is still un¬ 
clear, it should be kept in mind that DM scattering off stars 
will generally be much less efficient than in supergiant galaxies, 
as there are far fewer stars, if any, in the cores of the ultra faint 
dSphs that are heavily DM-dominated. Indeed there is increas¬ 
ing observational evidence that dwarfs contain central massive 
black holes. The AGN fraction may be as large as ~ 10% for the 
lowest mass nearby infrared-detected dwarfs and ~ 1 % if opti¬ 
cally selected with active AGN found for black holes in the 
mass range 10^^ lO^M© 31 . X-ray selected dwarfs are rarer, 
with central massive BH candidates amounting to ~ 0.1% [3^ . 
Some degree of theoretical plausibility is added, in addition to 
the hierarchical structure arguments cited previously, by the 
suggestion that IMBHs in dwarfs are a crucial element of early 
feedback, in terms of ejecting baryons, and reducing the baryon 
content of the host galaxy [^, an outstanding problem for 
MWG-like galaxies where the missing baryons are presumed 
to have been heated and ejected beyond the virial radius by 
unknown processes 34 . Gonstraints on minispikes have been 
studied in the case of WIMPs primarily producing a secondary 
gamma ray flux [3^ 36 , however dSphs in general have been 


poorly exploited for DM searches by neutrino detectors |37]. 


Remarkably, indirect detection of decaying DM via neu¬ 
trino detection becomes compatible with gamma ray signals 
for DM masses above ~ 1 TeV [^, and might be relevant for 
heavy DM candidates as in [^. One can however expect that 
future improvements to neutrino telescopes, for example with 
more closely spaced strings similarly to an enlarged DeepGore- 


like detector [^, will improve sensitivity considerably down 
to ~ 100 GeV. Here, we explore a SUSY thermal DM scenario 
that predicts enhanced branching ratios to monochromatic ly 
signals at these energies and illustrate the detection capabili¬ 
ties of future neutrino telescopes. 

The focus of this paper is to analyse neutrino lines origi¬ 
nating from annihilating DM. As a concrete example of DM 
annihilation to neutrino final states, we take a supersymmetric 
scenario. The model consists of the Minimal Supersymmet¬ 
ric Model (MSSM) augmented by a RH neutrino superfield, 
which provides a mostly right-handed (RH) mixed sneutrino 
as the lightest supersymmetric particle (ESP) [^[^. This 
model is well motivated as it addresses two basic problems: 
the origin of neutrino masses and the nature of DM. Purely 
left-handed (LH) sneutrinos as given by the MSSM have been 
studied in EDI as DM candidates, however due to their large 
coupling to the Z boson, they are excluded from direct detec¬ 
tion searches as a dominant component of DM. By introducing 
RH chiral superfields, the LH sneutrinos can mix with its RH 
counterparts and lead to the required relic density while 
fulfilling the direct search exclusion bounds. On the other hand, 
the RH fermionic field, a RH neutrino, provides the mechanism 
to generate neutrino masses and a connection between DM and 
neutrino phenomenology. 

Several possibilities for neutrino mass mechanisms and their 
consequences for sneutrino DM phenomenology have been in- 
vestig ated, fo r instance Dirac masses [4Q|44 4^ , seesaw mecha¬ 
nisms |44|47| p^ and inverse seesaw mechanisms |50|51| . In this 
paper, we consider the simplest possibility of Dirac masses, e.g. 
without lepton-violating terms, called MSSM+RN as in 
|53| . Neutrino lines produced by various neutrino mass gen¬ 
eration mechanisms have been analysed in in an effective 
field theory approach without however enforcing the DM con¬ 
straints. 


In our case, the sneutrino is the dominant DM compo¬ 
nent and compatible with the exclusion bound by LUX 
for direct searches. We show that sneutrino pair annihilations 
produces a sharp and enhanced monochromatic neutrino line 
at tree level, besides the usual secondary neutrino flux. As 
will be explained below, apart from the mass and mixing an¬ 
gle of the lightest sneutrino, the monochromatic neutrino line 
is mainly sensitive to the mass spectrum and nature of the 
neutralino. Thus, it provides an excellent complementarity to 
collider searches, in case the rest of the MSSM+RN spec¬ 
trum is beyond the reach of the LHG. Other indirect detec¬ 
tion prospects for sneutrinos have been studied in [S 46 
while Refs. |S 45 56 57 present a detailed analysis of the 
neutrino flux or line from the Sun or from the Earth. In gen¬ 
eral, the MSSM+RN LHG signatures can be quite distinct 
from those of the conventional neutralino LSP, as has been 
^ 58 ^ and in in the framework of sim- 


discussed in 

plified model spectra (SMS). Erom this latter analysis, which 
takes into account an extensive sampling of the MSSM+RN 
parameter space, it has emerged that there are almost always 
parameter combinations such that the limits from SMS can be 
avoided. It was shown that the direct detection experiments 
provide an excellent complementarity to the LHG searches for 
thermal sneutrino DM |53|. 


At a practical level, we consider a subset of points anal¬ 
ysed in [^. The MSSM+RN have 13 free parameters (gaug- 
ino masses, scalar masses for LH sleptons and squarks, RH 
sneutrinos, RH charged slept ons and squarks, trilinear con- 
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plings for all the scalar sector and finally the ordinary param¬ 
eters describing the Higgs sector) whose input value is given 
at the Grand Unification scale ~ 10^® GeV. The MSSM+RN 
SUSY particle spectrum is computed with the code SoftSusy, 
appropriately modified to adapt to micrOMEGAS_3.6 and 
micr0MEGAS_4.1 for relic density and particle differential 
spectra computations. We use a nested sampling algorithm 
MultiNest_v3 .2 [65[|66| to sample efficiently these free pa¬ 
rameters. In the model likelihood calculations, we include all 
DM constraints: the sneutrino satisfies relic density constraints 
from Planck measurements and the elastic scattering cross- 
section off nuclei is compatible with the 90% confidence level 
(GL) of the LUX exclusion bound [^. All details about the 
sampling procedure, the constraints and measurements imple¬ 
mented in the likelihood function, and the LHG phenomenol¬ 
ogy are provided in [^ . 

The rest of the paper is organised as follows. In the next 
section we discuss the modelling of the DM density spikes in 
SMBHs, the dwarf spheroidal galaxies and the DM minispike. 
In Sec. we study a simplified model for the monochromatic 
neutrino line to extract the relevant dependence on the the¬ 
oretical parameters. Section describes briefly how neutrino 
flavor oscillations modify the line signal. More importantly, it 
states the setup for an ideal neutrino telescope capable of de¬ 
tecting the monoenergetic neutrino lines in the 0(100) GeV 
energy range. In Sec. we provide details of the MSSM+RN 
model in which the sneutrino is a successful DM candidate and 
discuss the expected flux of monochromatic neutrino lines. We 
give our conclusions in Sec. 


2 Neutrino flux from WIMP annihilation in 
astrophysical sources with high DM density 


Sneutrinos as DM candidates can annihilate with a small but 
non negligible probability into the Standard Model (SM) par¬ 
ticles producing several possible final states, along with neutri¬ 
nos. The expected neutrino differential flux in an astrophysical 
source is given by: 


_ 1 ^2 {crv) dNjy ^ 


( 1 ) 


The first part of Eq. depends only on the WIMP model. In 
particular ^ (relevant for subdominant dark 

matter), (av) is the total annihilation cross-section of sneutrino 
pairs at the present epoch, ?n|^ is the sneutrino mass and 
dNjy/dE the differential neutrino spectrum per annihilation 
event, as a function of energy E: 


dJV^ ^ iBl 5{E - me) line, 

dE secondary. 


2.1 Dark matter spike in supermassive black holes 


We begin with supermasssive black holes in massive galaxies. 
The MWG SMBH is not the optimal target, mostly because 
the presence of a nuclear star cluster guarantees that scatter¬ 
ing will soften any DM spike, the scattering time-scale being 
of order of a Gyr [^. Of course continuing growth of the BH 
may partially compensate this, but the detailed evolution is 
complicated by the competition between growth of the nuclear 
star cluster by both gas accretion and in situ star formation as 
well as infall of globular star clusters . A much better target 
is a SMBH in a massive spheroidal galaxy such as M87 or Gen- 
taurus A (GenA). In the former case the stellar heating time 
of any putative DM spike is of order 10^^ Gyr, and therefore 
presents the possibility of an excellent DM signal amplifier [26| . 

We then consider three optimised galaxy targets, known to 
host a SMBH in their core and we assume that a DM spike 
has formed and survived all possible disruption processes. For 
instance, we analyse M87, a supergiant elliptical galaxy in the 
constellation Virgo at a distance of 16 Mpe from the Sun. In 
M87, it is plausible that the scattering off stars has been ineffi¬ 
cient in erasing the spike because it is dynamically young. The 
relaxation time for M87 is estimated to be E ~ (!)(10^) Gyr due 
to the strong dependence on the velocity dispersion of the stars 
and the relatively low core density^ Hence, a spike formed at 
early times is much more likely to have survived galaxy dy¬ 
namics up to the present epoch in M87 than in the Milky Way, 
which has a relaxation time of ~ 2.5 Gyr. Similar reasoning 
holds for NGG1277 and GenA. We consider NGG1277, which 
is a lenticular galaxy in the constellation Perseus, as it may 
host one of the largest SMBH ever measured and is at 

a distance similar to M87 with respect to us. Finally, GenA is 
in the constellation Gentaurus and despite of hosting a smaller 
SMBH compared to M87 and NGG1277, it has the advantage 
of being one of the closest radio galaxies to the Sun. All rel¬ 
evant details about these three galaxies and their SMBHs are 
provided in Tab. 

The dark matter density spike for a SMBH consists of two 
components, ppi the DM density of the plateau and psp the 
spike profile. With the assumption of adiabatic growth of the 
DM spike, its density is given by p^ : 


f 0 


pm{r) = < 


Psp(r)Ppi 

Psp(r)+Ppi 

PO 




r < 4Rs 
4Rs <r < Rsp 

^ ^ Rsd : 


( 3 ) 


with Rsp being the spike radius. This expression relies on the 
assumption that the spike has grown from a DM density pro¬ 
file oc pQ{r/rQ)~^. The plateau density depends on the DM 
annihilation rate as: 


_ rup 

(crr’)tBH 


( 4 ) 


Bi is the fractional contribution into neutrinos of each channel 
k and diW. /dE is the number of neutrinos produced in the jth 
channel, [j For processes dominated by s-wave, see Sec. the 
particle physics part of Eq. ^can be computed independently 
of the astrophysical factor, ^Astro, which we discuss below. 


with ten the age of the BH with mass Mbh- The spike profile 
is 

Psp{r) = png-<{r) (^) > ( 5 ) 


^ The sneutrinos are not self-conjugate particles, hence the 
secondary neutrino differential spectrum in Eq. ^ should be 
rescaled by an additional factor 1/2. 


^ Dynamical heating by stars is inefficient when the dynami¬ 
cal relaxation time E in the BH core is larger than the Hubble 
time (~ 10^° yr). 
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Table 1. Parameters relevant for the analysis for each SMBH we consider. In all cases we assume Ibh — 10^° years, while black 
hole mass Mbh and distance D is taken from [^[^ . 


SMBH 

Mbh [Mq] 

Rs [pc] 

D [Mpc] 

Po [GeV cm 

^Astro [GeV" cm-"] 

Declination 

M87 

6.4 xl0“ 

6.1 X 10“^ 

16.4 

2.3 

3.5 X 10“ 

+ 12 ° 

GenA 

5.5 xlO^ 

5.3 X 10“® 

3 

9 xlO^ 

4.3 X 10"“ 

hL 

CO 

o 

NGG1277 

1.7 xl0^° 

1.6 X 10“® 

20 

495 

2.5 X 10^" 

+41° 


with 


2.2 Dwarf spheroidal galaxies 


o 

II 

ro J ’ 

(6) 

9j{r) « (l 

4Rsy 

r J ’ 

(7) 

Rsp — O!.jT0 

/Mbh ^ 

V POT’g J 

(8) 


Let us now turn our attention to ^Astro from dSphs. The ex¬ 
pected flux luminosity from dSphs is given by the so-called 
J-factor [T^ : 

$Astio = J= [ dn' f PdwB.vt{r{s,0))ds , (11) 

J J los 


and 7 sp = For all numerical evaluations we set ro = 

20 kpc similar to the Milky Way, — 0.1 and 7 = 1 (which 
correspond to a NFW density profile and gives a spike slope 
7sp = 7/3). The normalization po is fixed according to the pre¬ 
scriptions in and is in agreement with the values reported 
in [^ . 

Considering the galaxies hosting the SMBH as point sources, 
the astrophysical flux is given by: 

1 

^Astro — Aesc / dr 7* Pbh(^) ? (9) 

^ Jarq 

where D and Rs are respectively the distance to the Earth and 
the Schwarzschild radius of the SMBH. iVesc parametrises the 
fraction of neutrinos produced by WIMP annihilations that 
can actually escape from the BH. Numerical simulations tend 
to indicate that Aesc is always larger than 90% even close to 
the event horizon for a Kerr BH [71| , a value we take as a 
reference for the numerical analysis. 

The behavior of pBii{r) is shown in Fig. [T top panel, for 
a DM mass of 100 GeV and two different values of the to¬ 
tal thermally-averaged annihilation cross-section. For a cross- 
section close to the value of a standard thermal relic (i.e. 
~ 10 “^® cm^/s), the plateau region is clearly visible, while if 
(av) < 10“^° cm^/s the DM interaction is not strong enough 
to flatten out the spike for NGC1277 and M87. This means 
that below a certain value of {(Jv), it is possible to probe the 
entire DM spike leading to a maximal enhancement in the neu¬ 
trino flux, independent of the particle physics. This is not true 
for CenA, where the DM density is larger. In this case, even 
small cross-section values of 10~^° cm^/s flatten out the DM 
spike. We can approximate the DM density profile as: 


pBH{r) = 


4 


Ppi if Psp (t*) > > Ppi, 
Psp{r) if psp{r) << Ppi. 


( 10 ) 


In the former case, the astrophysical flux is proportional to 
the plateau density times R^/D^, while in the second case the 
dependence is ^Astro ~ plrl{0{Rs) + 0{Rsp)). Hence, in the 
first case, the luminosity gets enhanced for very massive even 
though distant SMBHs, while in the latter case this is not 
true anymore. This will be a relevant issue for our study, as 
discussed in Sec. [H 


with s being the distance along the line of sight and 0 the 
opening angle of the cone. For the J-factor, we use the values 
computed in for integrating over a 1° angle. This allows 

a better assessment of the various sources of uncertainties that 
affect the J-factorI3 

Numerical simulations suggest that IMBHs with masses 
between 10^ Mq — lO^M© may have grown and formed a spike 
in up to 10% of the cases for Draco-like dSphs [^, and this will 
be our underlying working hypothesis. To model this minispike 
hosted at the center of the dSph, we use a similar prescription 
for pdwarf(^) as for the SMBHs except for a modification in 
Eq-i po is now defined as the normalization for NFW density 
profile and ro is the scaling radius as usual. The modeling of 
Psp{r) is given by Eq. ^ with Rsp and pR obtained with the 
prescriptions in [^[^ . 

While about 25 nearby dSphs were known as of early 2015, 
recently imaging data from the Dark Energy Survey has led to 
the discovery of nine new Milky Way satellites in the Southern 
sky increasing the total number to 34. Among all these 

dSphs, we consider the 8 galactic dSPhs that have the best 
measured stellar kinematics, and 14 ultra faint galactic dSPhs, 
for which the kinematical data are more uncertain but however 
are interesting due to their proximity and sky coverage. The 
relevant numbers are reported in Tab.j^and the separation into 
Northern and Southern hemispheres will be useful for assessing 
the sensitivity of neutrino telescopes. 

Table lists all relevant parameters for the two dSphs we 
consider in detail, Draco in the Northern hemisphere and Retic¬ 
ulum H in the Southern sky, for two different IMBH masses, 
lO^M© and lO^M©. The total luminosity is given by the flux 
along the line of sight from the minispike plus the standard 
J-factor, Eq. EH The choice of these two dSPhs is dictated 
by the robustness of the J-factor against uncertainties ME- 
Reticulum H, in addition to being robust against uncertain¬ 
ties has one the largest J-factors among all dSPhs listed in 
Tab.|^ In the Northern sky, the largest J-factors are provided 
by Ursa Major H and Willman I; however these ultra faint 
galaxies possess either a velocity held indicating strong tidal 
disruption effects from the Milky Way or nonequilibrium kine¬ 
matics. Given its proximity and large DM content. Reticulum 


^ Uncertainties include for instance the choice of the DM 
density profile, the sample size and the quality of kinematic 


data; for a thorough discussion we refer to 73 
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Table 3. Minispike parameters, po and scaling radius f or Dra co and Reticulum II for two assumed BH masses, ro and po for 
Draco are taken from 


36 


while for Reticulumll we use 74 76]. In all cases we assume tsH = 10^^ years. 



dSph 

ro [kpc] 

po [GeV cm 

Rsp [pc] 

PR [GeV cm“'^] 

Mbh = 10‘‘(10^)M© 

Draco 

2.09 

0.99 

1.5 (0.15) 

1.3 X 10''(1.3 X lO"") 

Rs = 9.57 X 10“®(xl0"^^) pc 

Reticulum H 

4.28 

2.81 

0.63 (0.063) 

1.9 X lO'^O.g X 10®) 


Table 2. Distance and J-factors (1°) for each dSPh we consider 
in the analysis, from lUlZll 


dSph 

D [kpc] 

1 J(l°) [GeV^ cm-^] 

Northern sky 

Draco 

80 

2.11 X 10^“ 

Ursa Minor 

66 

1.24 X 10^® 

Sextans 

86 

8.09 X 10^^ 

Leo I 

250 

8.87 X 10^^ 

Leo H 

205 

1.37 X 10^® 

Northern sky (ultra faint) 

Segue I 

23 

2.06 X 10^'' 

Ursa Major H 

30 

1.87 X 10^° 

Segue H 

35 

1.72 X 10^® 

Willman I 

38 

4.75 X 10^® 

Goma 

44 

8.32 X 10^® 

Bootes I 

66 

6.07 X 10^® 

Ursa major I 

97 

6.79 X 10^® 

Hercules 

132 

1.99 X 10^® 

Ganis Venatici H 

160 

4.13 X 10^^ 

Ganis Venatici I 

218 

4.50 X 10^® 

Leo V 

180 

1.88 X 10^® 

LeoT 

407 

4.80 X 10^^ 

Southern sky 

Garina 

101 

1.05 X 10^® 

Fornax 

138 

7.07 X 10^^ 

Sculptor 

79 

4.30 X 10^® 

Southern sky (ultra faint) 

Leo IV 

160 

2.14 X 10^® 

Reticulum H 

30 

5.88 X 10^° 


II seems to be one of the most promising objects for WIMP 
annihilation searches. 

Figure shows the behavior of the DM density profile as 
soon as the DM spike sets in for Draco, in the bottom panel. 
The DM density is enhanced by a few orders of magnitude with 
respect to the case of the SMBHs (top panel). As opposed 
to the case of the SMBHs, small annihilation cross-sections 
cm^ s“^ are not able to probe the full spike, as clearly 
indicated by the presence of the plateau due to WIMP anni¬ 
hilation. The change in the BH mass in dSphs does not signif¬ 
icantly modify the picture: a change from Mbh = lO^M© to 
Mbh = lO^M© changes the cutoff at ARs by a factor of 10^. 
The overall larger DM density and the vicinity of the dSphs 
with DM minispikes partially compensates the smallness of Rs 
in Castro ~ Thus we expect to have a larger neutrino 

flux from the dSph + minispike with respect to the SMBHs. 

It is clear that observations of dSphs with minispikes pro¬ 
vide a more promising way to observe monochromatic neutrino 
lines than SMBH. There is of course an additional layer of un¬ 
certainty: SMBHS are known to be present in massive galax- 




Fig. 1. Top: The DM density profile distribution for the three 
SMHBs is shown as a function of the distance from the center, 
as labelled. Bottom: Same as top panel but for a BH mass of 
lO^M© in the Draco dSph. The NFW density profile is denoted 
by the solid gray line for reference. In both panels we show the 
density plateau for different values for the thermal averaged 
cross-section as labelled, for a fixed DM mass of 100 GeV. 


ies, whereas the occurrence of IMBHs in dwarfs or ultra faint 
dwarfs is dependent on theoretically-motivated but plausible 
speculation, with some modest observational support for the 
presence of IMBHs in the more massive dwarfs. Having com¬ 
pleted our discussion about the astrophysical part of Eq. 
we now turn our attention to the particle physics model under 
consideration and sketch the conditions necessary to get the 
monochromatic neutrino lines from WIMP annihilation. 
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3 Analysis of the monochromatic neutrino 
line 

3.1 MSSM+RN model 


The model we consider is the MSSM with the inclusion of the 
RH neutrino superfield (MSSM+RN) as defined in [^[40l[44] . 
We briefly illustrate the main features of the model for clarity. 
The superpotential for Dirac RH superfields is given by 

W = CijijiHiHf - + Vy Hi L]N^) , (12) 

where is a real and diagonal matrix in flavor space, from 
which the neutrino mass is = VuYJ^ . It should be noted 
that lepton-number violating terms are absent in this frame¬ 
work and their inclusion is left for future work. The RH scalar 
fields add new terms in the soft-breaking potential 

Hsoft = {MlY^ LYLi + - 

[eij{AYHtyy + + h.c.], (13) 

where the matrices M% and A^ are real and diagonal, M% — 
diag(m^fc) and A^/ — diag(Ai>^) (k — are the flavor 

index as in Sec.|^. The sneutrino mass potential is defined in 
the sneutrino interaction basis, = (/>£, N*), as 


v'‘ = + 

^ mass ^ 




(14) 


with the squared-mass matrix 


Mlr — 

my + eos(2/3) + ml ^ 


(15) 


+2 


Here, m^ = are the soft mass terms for the three 
SU(2) leptonic doublets supposed to be common to all flavors, 
tan/3 = Vujvd and — (246GeV)^, where Vu,d are 

the Higgs vacuum expectation values (vevs). 

The off-diagonal term determines the mixing between LH 
and RH fields. The Dirac neutrino mass mo can be safely ne¬ 
glected as it is very small. All the mixing is hence provided 
by the trilinear term. If it is aligned to the neutrino Yukawa 
(Af>^ = rjYjy), this term is negligible compared to the diagonal 
entries and no efficient mixing is provided. However, can 
in general be a free parameter 39 40 and can induce a siz¬ 


able mixing among the sneutrino interaction eigenstates. The 
sneutrino mass eigenstates are: 


i^ki 


- sin 0^ cos 0^ 
cos Or, sin Or, 
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iV'= 


(16) 


Notice that scalar trilinear term for the sneutrino is connected 
to the mixing angle as: 


sin2^i>. = V2 


Arj^ V sin /3 
(m| — ) 

^ ^k2 ^kl'' 


(17) 


The sneutrino coupling to the Z boson, which does not interact 
with SU(2 )l singlets, is largely reduced if the mixing between 
LH and RH is sizable. This is extremely relevant to have a 
good sneutrino DM candidate compatible with the LUX ex¬ 
clusion bound, see e.g., Refs. [^[^[^. The limits from 
direct detection experiments thus restrict the mixing angle for 
sneutrino DM. 


^0 

Xi 


Fig. 2. Feynman diagram leading to the monochromatic neu¬ 
trino line: a pair of sneutrinos is annihilating into Ur via t- 
channel neutralino (i — 1,...,4) exchange (there is an analo¬ 
gous contribution from pair of anti-sneutrinos annihilating into 
anti-neutrinos). 


The renormalization group equations are modified by the 


new singlet superfields N as described in 45|[59 . In particular 


the new trilinear term induces an additional running to the 
LH part, while this is the only correction for the RH soft mass. 
Assuming common scalar masses and trilinear couplings for all 
flavors and neglecting all lepton Yukawas but Y- in the RGBs, 
the sneutrino tan, mo ^^, ends up to be the lightest one among 
the three sneutrino flavors and hence the LSP, while muon and 
electron sneutrino remain mass degenerate. 

It should be noted that the mass splitting between mo^^ 
and is frequently smaller than 5 GeV, resulting in prac¬ 

tically degenerate sneutrino spectrum at the LHG. However on 
cosmological time scale the heavier sneutrinos will eventually 
decay into , hence they do not play a role at the present 
epoch for annihilation in astrophysical objects. The rest of the 
SUSY particles are unaffected by the presence of the RH scalar 
fields and their behavior is the same as in the MSSM. 


3.2 Simplified model for the neutrino monochromatic 
line 

For WIMP annihilation at the present epoch there are two 
possibilities for producing a neutrino flux. The neutrino flux 
is usually considered to be given by secondary neutrinos re¬ 
sulting from the annihilation of sneutrinos into SM particles 
{i.e. W^W~, ZZ,hh,tt etc.), which subsequently 

fragment and decay, leading to neutrinos. A second possibil¬ 
ity is direct tree level (whenever possible) annihilation into 
monochromatic neutrinos, which will be the main part of this 
paper. Figure [^illustrates the tree-level Feynman diagram rel¬ 
evant for sneutrino pair annihilation into a monoenergetic neu¬ 
trino line. It only involves sneutrino and neutralino parameters, 
hence to study the behavior of the line with respect to the the¬ 
oretical parameters of the model, we take a simplified model 
approach involving sneutrino LSP and the lightest neutralino 
(z = I) for simplicity (this will turn out to be a good assump¬ 
tion when studying the realistic MSSM+RN). We neglect the 
DM constraints for now. The features of neutrino lines as well 
continuum for a viable sneutrino DM will be discussed in Sec.|5] 

On top of the SM particle content, the simplified model we 
consider is constituted by the sneutrino LSP z>ri, its mixing 
angle 0 ^^, the lightest neutralino y? together with its compo¬ 
sition: the Bino and Wino fraction {Nn and A '12 respectively) 
and the Higgsino fraction (Yi 3 ,Yi 4 ). In total there are two 
free parameters in the sneutrino sector, the LSP mass mj^^^ 
and sin Or ^, plus 4 free parameters from the electrowikino sec¬ 
tor. 

The process responsible for the monochromatic neutrino 
line is f'nf'ri ^ via t-channel exchange of a neutralino 
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Fig. 3. Left: The thermal averaged cross-section as a function of the sneutrino mass for each possible sneutrino annihilation 
channel. The sneutrino is assumed to be pure LH. For fixed mass of neutralino and mixing angles in the neutralino and sneutrino 
sectors, we vary the mass of the sneutrino. Right: The thermal averaged cross-section as a function of the sneutrino mixing 
angle; the trilinear coupling is fixed at 60 GeV and the SUSY masses are fixed as labelled. For fixed neutralino and sneutrino 
masses and mixing in the neutralino sector, we vary the mixing angle of the sneutrino. 


(see Fig.|^. The dominant term for the thermal averaged cross- 
section is s-wave and given by, see also [^: 




Ck 


+ Ck 


Xi 


Stt 


(m|o +m?J 


(18) 


with 


Cp^ = --sm0iy^YjyNi4 , 

1 a ^ AT ^ e COS OwNil sin , 

Cpl = --cos0u^Y^Ni4 H - ^ ^ , + 

2 2V2{1 - sin^ Ow) 


(19) 


e sin^p^A^i 2 sinOo 


eNi 2 sin 0^ 


2\/2{l — sin"^ Ow) 2\/2 sin Ow — sin"^ Ow) 

Yjy is the neutrino Yukawa coupling, e the electric charge 
and sin Ow the Weinberg angle. Since the exchanged particle 
is Majorana the annihilation process involves as well two anti- 
sneutrinos, which will give rise to two monochromatic anti¬ 
neutrinos. It is worth noting that there is a velocity dependent 
subdominant term in (crv), which can result in enhancement of 
the neutrino line in case of rotating BHs as mentioned earlier. 
We neglect this possibility for the current work. 

From Eq. 18 we notice that {av) is directly proportional to 
the mass difference between sneutrino and neutralino. Indeed 
for m^o » rriu^^ the thermally-averaged cross-section simpli¬ 
fies to oc l/m?o- From Eq. 

Xi 


19 


it is clear that the 


Higgsino component does not contribute to this process as it 
suppressed by the neutrino Yukawa coupling. This also implies 
Cph ~ oQln order to illustrate the dependence of annihilation 
cross-section and the differential spectrum of the neutrino line, 
we vary the physical masses of the neutralino and sneutrino as 
well as their mixi ng a ngles, independently of each other. The 
dependence of Eq. h^on the Bino and Wino fractions of Xi is 


^ Cp^ ~ 0 is an artifact of this particular model of sneutrino 
LSP, with the Dirac mass term only. In the case of Majorana 
masses, one could get a non-trivial contribution from Cp^. 


displayed in Eig.[^ together with the other annihilation chan¬ 
nels of the sneutrino. The neutrino line final state is the only 
one that depends on other SUSY particles while all the other 
channels involve couplings and masses of SM particles or that 
of additional Higgs bosons, once the sneutrino mixing angle is 
fixed. In Eig.l^ left panel, we choose a purely LH sneutrino {i.e. 
|sin^i>^| = I) and plot (av) for < m^o. As soon as the 

sneutrino mass approaches the neutralino mass {av) decreases, 
while the overall normalization is dictated by the neutralino 
composition. Clearly the Wino case is enhanced with respect 
to the pure Bino case, as it receives contributions from both the 
SU(2 )l and U(I) couplings. As soon as kinematically allowed, 
the annihilation channels into W~W~^ and ZZ open up, me¬ 
diated by s-channel exchange of n Z boson for LH sneutrinos. 
The hh final state comes from the scalar quartic vertex with 
the sneutrinos, while for the whole mass range shown the anni¬ 
hilation channel into fermions is always open (we only show the 
bb contribution for simplicity). The bb final state receives large 
contributions from the Higgs pole and a small enhancement at 
the Z pole. We do not plot the continuum for Wino case as it is 
independent of the neutralino nature. The right panel exhibits 
the dependence of (av) on the sneutrino mixing angle. is 
directly proportional to sin , hence for purely RH sneutrino 
the annihilation into monochromatic lines is not possible. The 
relative weight between different annihilation channels strongly 
depends on the mixing angle, the contribution from Higgs s- 
channel becomes important relative to the Z boson s-channel 
exchange as far as sin decreases. The annihilation channels 
involving the Higgs boson are also proportional to the value of 
the scalar trilinear coupling, which appear in the Higgs cou¬ 
pling, see Eq. |I3| Those considerations find further support 
from Eig. where the differential neutrino spectrum per an¬ 
nihilation event is shown. The first panel illustrates the de¬ 
pendence of dN/ (red) and of the maximum of the Ur 
continuum (blue) as a function of the sneutrino mass, assum¬ 
ing again a LH LSP. The line is enhanced over the continuum 
over all of the mass range, except at the Higgs pole, when 
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Fig. 4. Left: Differential neutrino spectrum per annihilation event as a function of purely LH sneutrino mass. The red curves 
denote the differential neutrino line spectrum, while the blue curves stand for the maximum of the continuum energy spectrum. 
Both Bino and Wino neutralino cases are considered with solid and dashed lines respectively. Center: Differential neutrino 
spectrum per annihilation event as a function of the sneutrino mixing angle. The mass spectrum is fixed as labelled and the 
neutralino is pure Bino. Two different cases for the trilinear scalar coupling are considered with solid and dashed lines. Right: 
Same as the central panel for Wino like neutralino. In all the three panels we vary the parameter on x-axis while keeping the 
rest of the parameters fixed. The mixing angle is varied independent of the trilinear coupling. 


sneutrino annihilation into bb dominates the thermally aver¬ 
aged cross-section. The suppression of the differential neutrino 
line spectrum is negligible at the Z pole. The continuum seems 
to depend on the neutralino nature, however this is an artifact 
of the overlap of the Bino and Wino cases for the line. Since 
the Wino contribution is larger, then the neutrino continuum 
differential spectrum decreases. dN /for both Bino and 
Wino neutralino overlap until the opening of the W^W~ final 
state, after which the Bino case has a smaller contribution. The 
central and right panels show the dependence of dN/dEjy on 
the sneutrino mixing angle. The line intensity (red lines) drops 
abruptly for sin 0, while it reaches a plateau value for 

sizable mixings and is larger than the continuum. We show the 
behavior for two different values of the scalar trilinear term: as 
expected the larger the smaller the contribution of the line 
with respect to the continuum. This holds for both the Bino 
(central panel) and Wino (right panel) cases. 

In summary, the conditions to have a sizable contribution 
from the process f'nf'ri with respect to h'nh'ri anni¬ 

hilations are: 1) sizable mixing angle, 2) small 3) sizable 
mass gap between LSP and lightest neutralino and 4) relatively 
light mass spectra. 

In the case of the neutralino LSP, the neutrino lines are fur¬ 
ther suppressed because of the Major ana nature of the initial 
state. The process is analogous to Fig. with the neutralinos 
exchanging a sneutrino in t-channel, with an expression equiv¬ 
alent to Eq. up to smaller normalization factors. Notice 
however that in this case the sneutrino can be purely LH and 
the line is not bound to be a monochromatic Ur line but could 
involve all three neutrino flavors. Because of the smallness of 
the line in the case of the neutralino LSP, we do not investigate 
any further the possibility (see Ref. for details). 


4 Detection of monochromatic neutrino lines 

The neutrino spectra produced by DM annihilations at the 
source (SMBHs or dSphs) position are shown in Fig. left 
panel. The neutrino with r flavor has a sharp line at the DM 
mass and a negligible continuum spectrum (red), contrary to 
the case of Ue and (blue and green). The Ue and only 


contribute to secondary neutrino flux and do not produce line. 
These spectra however are modified by neutrino flavor oscilla¬ 
tions as follows. 

4.1 Neutrino flavor oscillations in short 

The propagation of neutrinos produced by sneutrino annihila¬ 
tions is affected by flavor oscillations, first in vacuum, when 
they travel from the dSphs or the SMBHs to the Earth, and 
eventually by oscillations in matter when crossing the Earth 
to reach the neutrino telescope. 

We follow for the implementation of flavor oscillations 
in vacuum in the long baseline limit, as the flavor oscillation 
length is negligible compared to the distance of the astrophysi- 
cal source. As for flavor oscillations in matter, we use the results 
from [8Q| together with their choice of oscillation parameters 
for coherence 0 

The neutrino spectra we use in the rest of the analysis 
assume normal hierarchy for neutrino masses and a maximal 
crossing of the Earth (zenith angle cosS = —1). This last 
condition is more suitable to produce up-going events into 
the neutrino detectors. An example of the neutrino energy 
distribution after oscillations is shown in the right panel of 
Fig. Ill with respect to the injection neutrino spectrum (left), 
the overall effect of flavor oscillations is the democratic re¬ 
distribution of the line in all three neutrino flavors, namely 
~ l/3diV/dS— (fc = e,Ai,T). 

We do not provide further details on this part of the discus¬ 
sion, as different assumptions on the neutrino mass hierarchy 
do not have significant impact on the results of our analysis. 
The results for anti-neutrinos, are similar to those for us. 

4.2 Future neutrino telescope set-up 

Small values for the mixing angle are closer to ‘realistic’ sneu¬ 
trino DM scenarios. It is indeed well known that pure LH sneu- 
trinos as DM are excluded by direct detection constraints, and 

^ These are up-to-date values and take into account a non¬ 
zero mixing angle between the neutrino mass eigenstates 1 and 
3, i.e. 0i3 = 0.15. 
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Fig. 5. Left: Differential neutrino spectra per annihilation event as a function of the neutrino energy, at the source. Right: Same 
as left but at detection, after including vacuum and matter oscillation thought the Earth; normal hierarchy and zenith angle 
cos (5 = — 1 are assumed. A very similar spectrum arises for the antineutrinos. The parameters for this point are rriu^^ =118 
GeV, m^o = 127 GeV, N12 = 0.996, iVn = 0.001, {av} = 2.15 x 10“^^cm'^s-S Bune = 0.38, Ap^ = 6.78 GeV, sin6»£>^ = -0.024. 


only very small mixing angle are allowed by LUX, as we will see 
in the next section. Already in the simplified model framework 
we can see from Fig. that small sinOu^ imply values of (av) 
below 10“^^cm^s“^. These values of {av) are orders of magni¬ 
tudes below current constraints arising from DM searches at 
the galactic center by neutrino telescopes. IceCube (IC) detec¬ 
tor holds the most stringent upper bound on the thermally 
averaged cross-section for DM masses between 10^ — 10^ GeV, 
to be less than {av) ~ 5 x 10“^^cm^s“^. 

In this section we propose a search strategy for monochro¬ 
matic neutrinos from DM annihilations and provide a guideline 
for the ideal set-up of future neutrino detectors. We focus on 
the energy range between 100 GeV to TeV DM mass scales, 
one of the most common ranges predicted by scenarios beyond 
the Standard Model, and assume that the values of {av) are 
below the reach of present DM detectors. Such small values of 
{av) are characteristic of the sneutrino scenario we use as a 
proxy for the monochromatic line search, but can be realised 
in several other WIMP scenarios as well. 

The IG neutrino telescope consists of 5600 digital optical 
modules (DOMs) arranged vertically along 86 strings embed¬ 
ded in a km^ of extremely transparent ice below the South Pole. 
The DeepGore array constitutes the central part of IG and con¬ 
sists of six strings arranged with closer horizontal spacing and 
instrumented with DOMs of higher sensitivity and closer ver¬ 
tical spacing along the string. DeepGore is able to reduce the 
energy threshold of IG from 100 GeV down to 10 GeV, however 
the energy resolution for these energies, without directional in¬ 
formation like in the case of searches for DM from the Sun, is 
estimated to be not more than 50%. In general, it should be 
noted that neutrino detectors such as IG and Antares , the 
largest neutrino detector in water in the Northern Sky, have 
indeed been designed and optimised for ultra high energy neu¬ 
trinos {Ejy > few TeVs) and hence are not sensitive to the 
range of masses considered in this paper. 

The detection of neutrinos occurs via the observation of 
the Gherenkov radiation emitted by the particles produced in 
the interactions between the neutrinos and the matter in and 
around the telescope. Particularly interesting are muons cre¬ 


ated by the charged current interactions of as they create 
long and relatively detectable tracks compared to other lep¬ 
tons. An even more clean signature is that of up-going muons, 
which provide a powerful discriminant to distinguish galac¬ 
tic and extra galactic neutrinos from atmospheric neutrinos. 
By considering these up-going events, current neutrino de¬ 
tectors have a rather high sensitivity to point sources in the 
opposite hemi^here emitting ultrahigh energy neutrinos. For 


instance, IG4^ M has a sensitivity down to E^jy ~ 2 x 
10 “^GeVcm“^^^Mo sources in the Northern sky, while Antares 


83 


can probe down to E 


in the Southern sky. 

We first assume that a future generation of neutrino tele¬ 
scopes can achieve similar sensitivities as IG40 and Antares to 
point sources, for neutrinos in the energy range considered in 
the paper. In other words neutrino telescopes with the same 
effective area as for > TeV for 100 GeV — TeV DM parti¬ 
cles would better perform for DM searchesj^This roughly im¬ 
plies detectors of the current size, however with an increased 
string configuration and granularity of the detector, similarly 
to DeepGore, to augment the energy resolution. Once these 
requirements are met, the most promising targets for WIMP 
searches are to look for DM in overdense regions, such as dSPhs 
and DM spikes from central BHs in galaxies. 

In our analysis, we provide the expected flux of monochro¬ 
matic muon neutrinos, to rely on up-going muon track detec¬ 
tion, from DM annihilations (notice that Ve and Vr will provide 
a similar flux of up-going events giving rise to shower events 
instead of muons tracks) and compare it with the neutrino 
telescope sensitivity. The IG40 ideal detector can thus probe 
M87, NGG1277 and the dSphs in the Northern sky, while the 
Antares-like detectors would probe GenA and the dSphs in the 
Southern sky. Very optimistically, we give the expected sensi¬ 
tivity of the upgraded version of IG, a 10 km^ detector called 


IG40 means that the analysis has been done with data 
taken by IceGube with 40 strings. 

^ The effective area is by definition the area for which detec¬ 
tion efficiency of a neutrino is 100 %. 












































10 


Chiara Arina, Suchita Kulkarni, Joseph Silk: Neutrino lines from sneutrino DM 


IceCube-Gen2 (IC-Gen2), after 5 years of running and for 


the water km^ detector KMSNet 85 in the Northern hemi 


sphere as upgrade of Ant ares. We assume that the energy res¬ 
olution is about 30%, as it is expected for 0(1) TeV neutrinos, 
and an angular resolution of 1°. 

Although our estimates of the telescope sensitivities are op¬ 
timistic, our aim here is to analyse the potential of future neu¬ 
trino telescopes for detecting neutrino lines. We do not embark 
on a detailed study of the detection prospects in this work, nei¬ 
ther do we attempt to give any details of the possible telescope 
design required in order detect such a signal. 


5 Sneutrino DM from MSSM+RN: 
realizations for the neutrino monochromatic 
lines 


In this section, we present the prospects for detection of mono¬ 
chromatic neutrino lines in the full MSSM+RN model. 

The input parameters for the MSSM+RN are non-universal 
gaugino masses, the two Higgs doublet masses, /x and B are 
respectively the mass term for the Higgs fields in the superpo¬ 
tential and soft breaking potential, Eqs. |12| and|13| the trilinear 
couplings for charged sleptons, sneutrinos and squarks and the 
soft scalar masses for LH slept ons and squarks, for RH charged 
sleptons, RH sneutrinos and squarks. All MSSM+RN input 
parameters are defined at the Grand Unification scale ~ 10^® 
GeV. These 13 free parameters are efficiently sampled with 
the nested sampling package MultiNest_v3.2 [^[^, based 
on Bayes theorem. As a result, the sampled points are dis¬ 
tributed accordingly to the posterior probability density func¬ 
tion and satisfy the DM constraints: relic density constraints 
from Planck measurements [43| and the elastic spin indepen¬ 
dent (SI) scattering cross-section is compatible with the 
90% GL of the LUX exclusion bound. All details about the 
sampling procedure, the constraints and measurements imple¬ 
mented in the likelihood function, and the LHG phenomenol¬ 
ogy are provided in [^. At practical level the relic density 
and direct detection are computed with micrOMEGAS_3.6 [^ , 
while we use micr0MEGAS_4. 1 [64| for the neutrino spectrum 
and indirect detection annihilation cross-section. 

In this section we study the sampled points which have 
small Ajy^ terms and reasonably large sneutrino mixing angles, 
both n ecessa ry conditions to have neutrino lines, as discussed 
in Sec. 3.2P Of course, all these points satisfy the constraints 
listed in |53| . 

The differential spectrum for per annihilation event is 
shown in Fig. in the sneutrino-neutralino mass plane (top 
panel) and as a function of the sneutrino mass and thermally- 
averaged annihilation cross-section at the present epoch (bot¬ 
tom panel). Gonsidering first the bottom panel it should be 
noted that the typical value of (av) is smaller than the usual 
value for the standard WIMP relic. Indeed the largest values 
start around I0“^^cm^s“^ and go down even below 10“^®cm^s“^ 
For the sneutrino mass range covered by the sample, these 
values of (av) would be too low to achieve the correct relic 
density via its annihilation only (i.e. for these values the sneu¬ 
trino would cause overclosure). In the early Universe there is 


® This sample corresponds to logarithmic prior probability 
density functions on all parameters. 



Fig. 6. Top: The differential energy spectrum per annihilation 
event of Ur at the source position is shown as a function of both 
the sneutrino mass and the neutralino mass. Bottom: Same as 
top as a function of sneutrino mass and annihilation cross- 
section. The color code for dN /is given by the color 
bar on the right hand side. 


one primary coannihilation process at work: the lightest neu¬ 
tralino coannihilates efficiently with the lightest chargino and 
achieves the correct relic density (x? xf ^ qq\ ZW ^,and 
X?X? ^ VF+VP- and xtxl ^ W^W-, qq, ZZ,...), which is 
subsequently transmitted to the sneutrino sector. This argu¬ 
ment is further supported by the top panel, where m^o is 


shown as a function of the sneutrino mass. The sneutrino and 
lightest neutralino are almost degenerate, except at the Higgs 
pole and for few other points. The proximity in mass between 
the DM and Xi tends to cause a suppression in the intensity 
of the monochromatic Ur line, denoted by the color bar and 
the third dimension in both panels. Indeed the largest values 
for dN/dEj^^^^^ are achieved for sizable mass differences be¬ 
tween the lightest neutralino and the LSP (reddish points). 
In the few points above 100 GeV where the mass gap is siz¬ 
able, the sneutrino can thermalize due to coannihilation pro- 
. cesses with and/or with the charged sleptons, most 

usually the f^. On the contrary, even if at the Higgs pole, 
the mass gap between sneutrino and neutralino is large, the 
channels h ^ bb,r~^r~,... dominate and suppress 

all the other sneutrino interactions. These features were ex- 
pec ted f rom the discussion on the simplified model for the line. 
Sec. 3.2 The behavior of the line intensity dN/dEu^^^^ is further 
investigated in Fig. where the dependence on the sneutrino 
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Fig. 7. Upper left: Differential energy spectrum for the Vr line 
(blue) and for the continuum (green) are shown as a function 
of the sneutrino mass. Right: Same as left as a function of the 
sneutrino mixing angle, for the Ut line spectrum only. Lower 
left: Differential Ur line spectrum as a function of the Bino (red) 
and Wino (blue) fractions of the lightest neutralino. Lower 
right: Same as top right as a function of the scalar trilinear 
coupling. In all panels the differential spectrum is given at 
source. 


mixing angle and neutralino composition is shown. We show 
the dependence on the Wino (blue) and Bino (red) fraction 
of Xi ill fh® lower left panel. Most of the lightest neutrali- 
nos are Wino-like because of the coannihilation requirements 
in the early universe, i.e. in order to be within 10% in mass 
with the lightest chargino [^. When the neutralino is Bino- 
like, the relic density is set by cohannihilation with sleptons. 
In any case, the Higgsino components of the neutralino are 
always negligible, which positively impacts the process giving 
rise to the monochromatic neutrino line. The sneutrino mix¬ 
ing angle compatible with the LUX exclusion limit in direct 
DM searches are given in the top right panel. As one can see 
these values are pretty small and ranges from practically 0 (i.e. 
almost pure RH sneutrinos) to the maximal allowed value of 
\ sin 0i}^\ < 0.04. As discussed for a couple of benchmarks in 
Fig. dN /depends on sinOo^ and tends to reach a 
plateau value. The mixing angle value at which the plateau 
sets in depends strongly on : larger values of the trilinear 
scalar coupling tend to delay the onset of the plateau to larger 
mixing angles. The values of the trilinear coupling are given 
in the lower right panel: the largest value of dN /are 
indeed correlated with the smallest values of Hence the 
constraint by the LUX experiment may suppress the overall 
line intensity because it selects on average mixing angles much 
smaller than the value required to reach the plateau value for a 
definite SUSY spectrum configuration. SI elastic scattering of 
sneutrino off nuclei arises via the exchange of a Z and h boson 
on t-channel. The Z channel is independent of the sneutrino 



mp [GeV] 


Fig. 8. Top: DM density plateau, Eq. predicted by the 
MSSM+RN, as a function of the sneutrino mass and of the 
annihilation cross-section (third direction denoted by the color 
code and color bar on the right). Also shown with horizontal 
lines are the maximal values of the spike profiles, Eq.[^ for each 
of the SMBHs used in the analysis, as labelled. Bottom: Same 
as top for Reticulum II hosting a Mbh = lO^M© in the inner 
center. The maximal value of the density spike is indicated by 
the gray horizontal line. 


mass and leads to crxe ^ I0“^^cm^. Since the LUX bound, 
for the mass ranges 60 to 500 GeV, ranges from I0“^^cm^ to 
6 X 10“^^cm^, the Z exchange should be suppressed by almost 
6 order of magnitude, which explains the smallness of sin 0 ^^, 
while the Higgs exchange is controlled by . 

The other important contributions to the neutrino final 
state are the secondary neutrinos, leading a continuum in the 
flux. In the first panel of Eig.j^we illustrate the relative weight 
of the line spectrum with respect to the Ur differential sec¬ 
ondary spectrum (green dots) as a function of the sneutrino 
mass. The origin of the continuum flux is clear: at the Higgs 
pole it is due to sneutrino annihilations into fermions, while 
above the lU, Z and h thresholds it receives contributions from 
W~^W~, ZZ and hh final states respectively (the kinematic 
edges for the latter final state and for the W~^W~ threshold 
appear particularly clean in the plot). The line signal can be 
enhanced over the continuum by on average two orders of mag¬ 
nitude, the enhancement being bi^er for light sneutrinos and 
light but not degenerate Xi (Fig. as expected. 

The differential spectrum for the monochromatic line dis¬ 
cussed until now is given at the source position. As described 
































12 


Chiara Arina, Suchita Kulkarni, Joseph Silk: Neutrino lines from sneutrino DM 


in Sec. |4.1| the effect of neutrino flavor oscillations due to prop¬ 
agation from the distant galaxies to Earth is to reduce the line 
intensity by roughly 1/3 and to distribute it democratically 
among all flavors. This determines the dN /at detec¬ 
tion. 

To estimate the expected neutrino flux, Eq. it is neces¬ 
sary to estimate ^Astro, which in turns depends on the particle 
physics model in case of DM spikes, Eqs. and because of 
the annihilation plateau that can smooth out the DM spike. 
Here it should be made clear that the values of {av) that enter 
into computation of the DM density plateau, ppi oc l/{av), are 
those reported in Eigs. [^and|^ Eor sneutrino DM, it is crucial 
to differentiate between the value of {av) in the early uni¬ 
verse and at the present time. Those two quantities can differ 
by several orders of magnitude due to coannihilation effects, 
as already discussed above. Since after freeze-out, all SUSY 
particles decay into z>ri, at the time of the SMBH or dSPh 
formation, the DM is constituted solely of sneutrinos, annihi¬ 
lating with small (av). Eigurej^ illustrates the value range for 
the DM density plateau for the sampled points, as a function 
of the sneutrino mass and (cri’) (denoted by the third direc¬ 
tion and color bar on the right). The top panel refers to the 
SMBH case and shows the maximal value of the spike profile, 
~ 5 X lO^^GeVcm”^, for M87 and NGC1277 (horizontal vi- 
olet/green line). This value has to be compared with pp\ and 
for ppi > psp the enhancement from the DM spike is maxi¬ 
mal and the density plateau can not form anymore. This hap¬ 
pens for annihilation cross-sections below 10“^^cm^s“^. The 
spike profile of GenA (red horizontal line) has a larger maxi¬ 
mal value and the DM density plateau forms for cross-sections 
up to ~ 10“^^cm^s“^. The case of a dSph (Reticulum H) with 
a IMBH of mass Mbh = 10^M© is shown in the bottom panel: 
since dwarf galaxies are dominated by DM, the density spike 
is larger than the SMBH case (yellow horizontal line) and the 
density plateau always forms unless (av) < 10“^^cm^s“^. The 
value of the minispike profile, Eq. is not affected quantita¬ 
tively by the choice of the mass of the central IMBH. 

Having quantified all parameters that enter in the flux es¬ 
timate, Eq. we can now proceed and examine the detection 
prospects for monochromatic neutrino lines as a function of 
the sneutrino mass (which also indicates roughly the neutrino 
energy passing into the detector), as illustrated in Eig. Let 
us first consider the observation of dSPhs without the pres¬ 
ence of DM spikes (first two panels, top left). The dSphs in 
the Northern sky can produce up-going muon events in an 
IG-like detector located at the South Pole: we show the flux 
expected by Draco (cyan), Ursa Major H (blue) and the stack 
of all dSphs, as listed in Tab. The expected flux coming 
from the dSph stack could be in the reach of a detector with 
the same sensitivity as IG40 for point source, while the flux 
from Draco could only be measured by the IG-Gen2 like de¬ 
tector. The coannihilation mechanism acting at freeze-out has 
the drawback of suppressing the overall flux from sneutrino 
annihilations at the present epoch, as the viable sneutrino DM 
configurations exhibit such small values of {av ). The usual de¬ 
pendence of the flux on the DM mass as l/m^M ameliorates 
slightly the prospects for detection of 60 — 100 GeV sneutri¬ 
nos. By playing the same game with the dSPhs in the Southern 
hemisphere, which can give rise to up-going muons in a water 
detector such as KM3NeT, we see that the main contribution 
comes from Reticulum H (the flux from the stack of dSPhs co¬ 
incides with the flux produced by Reticulum H, gray points). 


Since there are only a few known dSphs in the Southern Sky, 
the flux is rather small and below the sensitivity of KM3NeT. 
It should be noted however that the 68% GL J-factor uncer¬ 
tainty ranges from 6.71 x 10^® GeV^ cm“^ to 1.14 x 10^^ GeV^ 


74 , which implies an uncertainty in the flux of roughly 


a factor 2. 

The inclusion of a minispike in Draco or Reticulum H ame¬ 
liorates the prospects for detection considerably (first two pan¬ 
els, bottom left). A minispike due to an IMBH of lO^M© (dark 
brown) or 10^MQ(light brown) increases the flux by 4 or two 
orders of magnitude respectively in both dSphs. Both scenar¬ 
ios can be probed with our ideal neutrino telescopes in the 
configuration of IG40 or Ant ares, while next generation de¬ 
tectors would investigate more deeply the MSSM+RN. Even 
though the minispike can enhance the predicted flux by 4 or¬ 
ders of magnitude, it still lies far away from the reach of the 
present IG detector in its real configuration, because the point 
source sensitivity for neutrinos below 1 TeV is of the order 
10“^GeVcm“^s“^. In the right bottom panel, we show 
the dependence of the flux on the thermally-averaged cross- 
section. In the absence of a minispike in the inner core of 
Draco, only cross-sections of the order 10“^®cm^s“^ can be 
reached with the future IG-Gen2, as the flux scales linearly 
with respect to {av ). The DM minispike introduces a nontriv¬ 
ial dependence on {av); the current IG40 sensitivity might be 
able to reach values down to 10“^^cm^s“^, that could not re¬ 
alistically be probed otherwise, for certain SUSY parameters. 

As expected, the prospect for detection of the flux coming 
from SMBHs (right panel on top, as labelled) is less promising 
than the dSph case. This is caused by two phenomena: 1) the 
SMBHs are much distant objects than the dSphs, and 2) the to¬ 
tal flux depends only on the DM density spike and not anymore 
on the density plateau, hence the flux cannot take advantage 
of the R%/D^ enhancement (i.e. the fact that the SMBH is 
actually very massive), as already anticipated in Sec. 
are not considering this case any further. 


2.1 


We 


5.1 Complementarity with other WIMP searches 


Gamma rays from dSphs - The Draco example: 

The sneutrino DM considered in this work appears to be quite 
elusive as far as the indirect DM searches are concerned. Sneu¬ 
trinos annihilating in dSPhs and SMBHs produce a secondary 
spectrum of gamma rays, mainly from the decay of W^W~ or 
bb and their subsequent fragmentation into pions. Even assum¬ 
ing a branching ratio 100% into W~^W~ or into bb, the current 
Eermi limits from dSPhs cannot constrain any of our sampled 
points, as the 95% GL upper limit in {av) versus DM mass lies 

for 


87 


in the ballpark from 3 x 10 ^'^cm'^s ^ to 10 “"cm"s 
100 — 500 GeV DM mass. Under the same assumptions for the 
branching ratios, the gamma ray flux produced by sneutrino 
DM is also compatible with the limits set by the DM spike in 


M87 [^. 

In the left panel of Fig. [Tol we illustrate the complementar- 
ity between future gamma ray and neutrino telescopes. In par¬ 
ticular we consider the Gherenkov Telescope Array (GTA), 
a ground based telescope array sensitive to high energy gamma 
rays. GTA is already in construction and may start its first 
scientific run around 2016 with a partial array. Its expected 
point source sensitivity is denoted by the blue dotted line. We 
show the expected flux for the neutrino line and for the gamma 
rays (orange points) from Draco, in the absence of a spike. 
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Fig. 9. Expected flux for the neutrino line as a function of the sneutrino mass, for all astrophysical environments considered in 
the analysis, as labelled in each panel. The flux is computed for and is given at the detector position. The color code for the 
sampled points is as follows. [Upper panels] Left: Cyan stands for Draco, blue for Ursa Major II and yellow denotes the stack 
of all dSphs in the Northern Sky. Center: Gray stands for Reticulum II (equivalent to the stack of dSphs in the Southern Sky). 
Right: Red denotes CenA, while light green and violet stand for M87 and NGC1277 respectively. [Lower panels] Left and eenter: 
Effect of the minispike for two different BH masses, as labelled. Right: Expected flux for the neutrino line as a function of the 
thermal annihilation cross-section in Draco, with(out) minispike in brown (cyan). In all panels we show the sensitivity of current 
and future neutrino detectors for point sources, assuming that it holds as well for neutrinos with energies ~ 100 ^ 500 GeV. 


The expected differential gamma ray flux is computed with 
micrOMEGAS_4.1 for all our sample points. We have integrated 
it in the energy range from 500 MeV up to DM mass and con¬ 
voluted it with a J-factor of 0.1° [^, which is the expected an¬ 
gular resolution of CTA. Clearly the intensity of the secondary 
gamma ray flux is similar to the monochromatic neutrino line 
flux. This is relevant in the sense that neutrino telescopes can 
be very competitive with next generation of gamma ray tele¬ 
scopes, if their energy resolution is improved down to 100 GeV 
energy scale. 

DM Direct detection: 

In the right panel of Eig. we illustrate the complementar¬ 
ity with direct detection searches. The values of the SI elas¬ 
tic cross-section for xenon detectors are plotted versus the ex¬ 
pected neutrino line flux. Independently of the astrophysical 
assumptions (with minispike, light brown, or without, cyan), 
the configurations that lead to the largest monochromatic neu¬ 
trino fluxes are in the reach of XENONIT after two years 
of running (vertical green line). XENONIT is the first ton-scale 


detector for DM direct searches currently under construction 
and is supposed to release the first data in 2017. Even though 
direct detection searches are insensitive to the presence of min¬ 
ispikes in dSphs, an interesting complementarity arises when 
looking at the brown sample, which includes a central BH, 
Mbh = lO^M©, in Draco. Indeed parameter configurations 
which are below the XENONIT sensitivity and close to the 
neutrino background (red vertical line) can be probed with 
the ideal IG-Gen2 or KM3Net detectors in five years time. 

The hypothesis of having sneutrinos as subdominant DM 
components might produce enhanced dN /with respect 
to values obtained in our sample. Indeed subdominant sneu¬ 
trinos can have large mixing angles, as the LUX exclusion 
bound gets weakened by the factor ^ / G/ipi^nck (where 

G/ipianck is the value measured by Planck for the relic density), 
i.e. <fcrxg. However, it should be noted that the enhancement 
in the intensity of the line is compensated by a reduction in 
the flux by so the full gain in visibility for the line is not 
straightforward. 
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Fig. 10. Left: Expected neutrino line flux compared to the expected gamma ray flux (orange) from sneutrino annihilation in 
Draco, as a function of the sneutrino mass. The point source sensitivity of CTA is shown together with the neutrino telescope 
sensitivities as labelled. Right: Predicted SI cross-section versus line flux. We show the complementarity of indirect and 
direct DM searches with and without minispike in Draco as labelled. The expected sensitivity of XENONIT and the neutrino 
background are shown by the vertical lines, as labelled. 


Complementarity with LHC searches: 

Another interesting complementarity arises when considering 
the BSM searches at the LHC. As pointed out in Sec. the 
neutrino lines considered in this work, concern only the sneu¬ 
trino and neutralino sectors. Eor the continuum, the addi¬ 
tional Higgs boson contributions matter. If the rest of the 
MSMM+RN spectrum is very heavy, the prospects of observ¬ 
ing signatures of MSSM+RN scenario at the LHC decrease 
rapidly. Indirect detection constraints at this point will be ex¬ 
tremely useful in order to constrain or even identify the sneu¬ 
trino DM scenario. 

6 Conclusions 

In this article, we have investigated the detection prospects for 
monochromatic neutrino lines. We take advantage of having a 
line signature pointing at the DM mass, and of having neutri¬ 
nos as final-state particles. Indeed, neutrinos are particularly 
clean signatures as they propagate freely and undisturbed from 
the source to the detector, contrary to gamma rays which might 
undergo absorption processes and may not leave the most cen¬ 
tral region of the galaxy e.g. close to the BH, especially in the 
case of massive galaxies. 

We have studied the detection prospects of such monochro¬ 
matic neutrino lines in two distinct astrophysical systems with 
a particularly high DM density: I. the dwarf spheroidal galaxies 
with and without a DM spike hypothesis, and 2. the DM spikes 
in SMBHs. These environments are amongst the best probes in 
order to search for DM and have already set very stringent con¬ 
straints in terms of gamma ray fluxes from DM annihilation. 
Remarkably, neutrino telescopes chafe the potential of doing at 
least as well for certain classes of supersymmetric DM models. 

As a particle physics candidate, we considered MSSM+RN, 
a well motivated and viable supersymmetric model with sneu¬ 
trino DM. The RH field is a necessary addition to the MSSM 
to provide both a successful DM candidate, a mostly RH sneu¬ 
trino, and to generate a Dirac neutrino mass term. The neu¬ 
trino lines in this model are a direct outcome of the annihilation 


of sneutrinos at tree level. Dark matter in dwarf galaxies spiked 
by IMBH provides an incredibly powerful means of probing 
low annihilation cross-sections well below 10~^®cm^s“^ that 
are otherwise inaccessible by any future direct detection or 
collider experiment. Eor instance, this signature depends only 
on a few SUSY parameters and in the case of heavy SUSY 
mass spectra, potentially out of reach of LHC, such neutrino 
lines can be a complementary signature to direct searches at 
the LHC. 

Sneutrino DM turns out to be very elusive in terms of in¬ 
direct detection. Sensitive neutrino detectors in ice or water 
are required to probe the monochromatic neutrino fluxes from 
dwarf spheroidal galaxies. The ideal set-up would be a detec¬ 
tor equivalent to IceCube but with increased granularity for 
the string and optical modules, in order to achieve the same 
resolution for (!I(300) GeV DM similar to those for ultra high 
energy neutrinos from point sources. Eor reasonable masses of 
the IMBH in the center of the dSphs, the expected flux of neu¬ 
trinos is similar to the gamma ray flux. In this case, neutrino 
detectors would be competitive with gamma ray telescopes, 
such as CTA, and add further complementarity to the multi¬ 
probe DM search approach. 

Erom a model building point of view there might be ways 
of increasing the tree level annihilation into neutrinos. Indeed 
for the Dirac neutrino mass term, the neutrino Yukawa cou¬ 
pling is very small and hence the line is sensitive to the LH 
component only, thus there is no tight connection between the 
neutrino line and neutrino mixing. The mixing angle in the 
sneutrino sector is instead limited by constraints from the LUX 
direct detection experiment. It may however be different, for 
instance, for the low-scale seesaw mechanism or seesaw type-I, 
where the neutrino masses are not generated by small neutrino 
Yukawa couplings. A more general investigation of such MSSM 
extensions and their impact on neutrino line detection might 
be interesting and is left for future work. 

The possibility of detecting monochromatic neutrino lines 
from DM annihilations in several nearby astrophysical systems 
provides an especially appealing case to look for a smoking gun 
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signature of a plausible dark matter candidate, the sneutrino, 
that is predicted by supersymmetry. The next five years of 
LHC running will be crucial for searching for supersymmetric 
signatures, and the complementary approach offered by future 
neutrino telescopes may be equally rewarding for exploring the 
sneutrino as a DM candidate, and provide an excellent oppor¬ 
tunity for fleshing out the multi-messenger study of the nature 
of DM. 
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